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Abstract: We quantify our current knowledge of the size and flavor structure of the matter 
effects in the evolution of atmospheric and long-baseline neutrinos based solely on the 
analysis of the corresponding neutrino data. To this aim we generalize the matter potential 
of the Standard Model by rescaling its strength, rotating it away from the ee sector, and 
rephasing it with respect to the vacuum term. This phenomenological parametrization 
can be easily translated in terms of non-standard neutrino interactions in matter. We 
show that in the most general case, the strength of the potential cannot be determined 
solely by atmospheric and long-baseline data. However its flavor composition is very much 
constrained and the present determination of the neutrino masses and mixing is robust 
under its presence. We also present an update of the constraints arising from this analysis 
in the particular case in which no potential is present in the e/i and er sectors. Finally we 
quantify to what degree in this scenario it is possible to alleviate the tension between the 
oscillation results for neutrinos and antineutrinos in the MINOS experiment and show the 
relevance of the high energy part of the spectrum measured at MINOS. 
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1. Introduction 

It is now an established fact that neutrinos are massive and leptonic flavors are not sym- 
metries of Nature [1,2]. In the last decade this picture has become fully proved thanks to 
the upcoming of a set of precise experiments. In particular, the results obtained with solar 
and atmospheric neutrinos have been confirmed in experiments using terrestrial beams of 
neutrinos produced in nuclear reactors and accelerators facilities [3]. The minimum joint 
description of all the neutrino data requires mixing among all the three known neutrinos 
(z^e, u^, Vt), which can be expressed as quantum superpositions of three massive states vi 
{i = 1, 2, 3) with masses mj. This implies the presence of a leptonic mixing matrix in the 
weak charged current interactions [4, 5] which can be parametrized as: 

C12C13 S12C13 Sl3e"*'^'^^\ 

CArac = ( -S12C23 ^ Ci2.Sl3S23e*''*^P C12C23 - Sl2Sl3-S23e*^'^^ C13S23 , (1-1) 

S12S23 - Cl2Sl3C23e*^^^'' -C12S23 - Sl2Sl3C23e*^CP C13C23 / 

where Cij = cos 6ij and Sij = sin 6ij. In addition to the Dirac-type phase 5cp, analogous to 
that of the quark sector, there are two physical phases associated to the Majorana character 
of neutrinos, which are not relevant for neutrino oscillations [6, 7] and which are therefore 
omitted in the neutrino oscillation analysis. 

In the simplest quantum-mechanical picture, flavor oscillations are generated by the 
kinematical Hamiltonian for this ensemble -ffyac which in the flavor basis {vg, ^t) reads 

Hy^c = Uy^cDy^Ul^ wlth Dy^ = diag(0, Amli, Am|i) (1.2) 

The quantities Am^i, \Am'^i\, 612, and ^23 are relatively well determined by the analysis 
of solar, atmospheric, reactor and long-baseline (LBL) experiments, while only an upper 
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bound is derived for the mixing angle ^13 and barely nothing is known on the CP phase 
6cp and on the sign of Am|^ [8-1 1] . Given the observed hierarchy between the solar and 
atmospheric mass-squared splittings there are two possible non-equivalent orderings for the 
mass eigenvalues, which are conventionally chosen as 



As it is customary we refer to the first option, Eq. (1.3), as the normal scheme, and to 
the second one, Eq. (1.4), as the inverted scheme; in this form they correspond to the two 
possible choices of the sign of Aml^.^ 

The flavor evolution of this neutrino ensemble is also affected by the difference in the 
matter potential induced by neutrino-matter interactions when it propagates in a back- 
ground of sufficiently dense matter through the so-called Mikheev-Smirnov-Wolfenstein 
(MSW) mechanism [12, 13]. Within the context of the Standard Model (SM) of par- 
ticle interactions, this effect is fully determined and leads to a matter potential which 
for neutral matter is proportional to the number density of electrons in the background 
Ne{r), V = \/2GFNe{r), and which only affects electron neutrinos. The evolution of 
the ensemble is then determined by the Hamiltonian = Hy^c + -f^mat' with iJ^at = 
^/2GFA^e(r)diag(l,0,0). The presence and magnitude of this potential in the propaga- 
tion of solar neutrinos can be tested in solar neutrino experiments (and in combination 
with KamLAND), and, as pointed out in Ref. [14], it agrees well with the SM prediction. 
Conversely the effects associated with non-standard forms of the matter potential in so- 
lar neutrino propagation have been studied and constrained by the analysis of solar and 
KamLAND data [15-26]. 

Matter background effects also affect the evolution of atmospheric and LBL neutrinos 
when traveling in the Earth. However standard matter effects for atmospheric and LBL 
neutrino oscillations arc suppressed when compared to solar neutrinos. This is so because 
atmospheric and LBL neutrinos are dominantly f^'s at the source, the Earth matter den- 
sity is smaller, and the characteristic length of the neutrino path in matter is shorter. 
Nevertheless the presence of non-standard matter effects are known to be relevant when 
discussing the precise determination of the oscillation parameters in the present and near 
future atmospheric and LBL experiments [27-52]. 

In this article we address our current knowledge of the size and flavor structure of 
the matter background effects in the evolution of atmospheric and LBL neutrinos based 
exclusively on the analysis the present data. In order to do so we introduce in Sec. 2 a 
generalized phenomenological parametrization of the matter potential allowing for rescaling 
of its strength from the SM prediction, rotation from the ee sector, and rephasing with 
respect to i?vac- We also discuss the connection with the matter potential induced by 
non-standard neutrino interactions (NSI) in matter which provide a well-known theoretical 
framework for this phenomenological parametrization. In Sec. 3 we present the results of 
our determination of this generalized matter potential from the analysis of atmospheric, 
LBL and CHOOZ experiments. In particular we will show that in the most general case, 
the strength of the potential cannot be determined by these data. However its flavor 
composition is very much constrained. We will also show that the present determination 

^In this convention the angles 9ij can be taken without loss of generality to lie in the first quadrant, 

dij G [0,7r/2], and the CP phase Sep € [0,2n]. 



Amli < (ATn|2 ^ Am|i) with (Am|i > 0) ; 
Am|i <C |Am|i ~ Am|2| with (Am|i < 0) . 



(1.3) 
(1.4) 
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of the neutrino masses and mixing is robust even in this generalized scenario. In Sec. 4 we 
will discuss the particular case in which no potential is present in the e/j, and er sector. 
In this case, a strong bound on the strength of the potential arises from the analysis of 
atmospheric and LBL neutrinos. We will update the results of Ref. [3, 36] and we will 
revisit the claims that within this scenario it is possible to alleviate the tension between 
the oscillation results for neutrinos and antineutrinos in the MINOS experiment. In Sec. 5 
we will draw our conclusions. Some technical details on the choice of parametrization of 
the matter potential are given in two Appendices. 

2. Formalism 

In the three-flavor oscillation picture, the neutrino (and antineutrino) evolution equation 
reads: 

















(2.1) 











where x is the coordinate along the neutrino trajectory and the Hamiltonian for neutrinos 
and antineutrinos is: 



= Hyac + -ffmat and = (i^vac " -f^mat)* > (2-2) 

with -ffvac given in Eq. (1.2). In the Standard Model .ffmat is fully determined both in 
its strength and flavor structure to be = \/2Gi;'A^e('") diag(l) 0, 0). In this work we 

generalize the form of the matter potential to be 

' Qrel = diag(e^"^e*"^e-^"l-^"2)^ 

C/mat = Rl2(.'Pl2)Rn{Vl3) , (2.3) 

D^^t = V2GFNe{r) diag(e, 0, 0) 

where we denote by Rij{ipij) a rotation of angle ipij in the ij plane. In this parametrization 
e represents a rescaling of the matter potential strength, 9912 allows for projection of the 
potential into the flavor and 1^13 allows for its projection into the Ut flavor. The two ad- 
ditional phases a\ and 02 included in Qrei are not a feature of neutrino-matter interactions, 
but rather a relative feature of the vacuum and matter term: they would become unphysi- 
cal if any of the two terms weren't there. As described in Appendix B one can always take 
< Oij < 7r/2 and Am|]^ > and the neutrino mass hierarchy is accounted by the sign of e: 
positive for normal ordering, negative for inverted ordering. For the case of real potential 
(ai = a2 = 0) one must consider —tt/2 < (pij < 7r/2, whereas for complex NSI (free ai 
and 0:2, < ttj < 27r) it is enough to assume < fij < 7r/2. In this parametrization the 
standard Sv oscillations in matter are recovered in the limit e = ±1 and (pi2 = (pi3 = 0, 
where e = +1 corresponds to the normal ordering and e = — 1 corresponds to the inverted 
ordering. 

As explained in Appendix A (see Eq. (A.l) and below) Eq. (2.3) is not the most general 
parametrization for the matter potential, since it has built-in the assumption that two of its 
eigenvalues are equal. In Ref. [37] it was shown that strong cancellations in the oscillation 
of atmospheric neutrinos occur in this case. Consequently despite our phenomenological 
parametrization in Eq. (2.3) represents only a subspace of the most general parameter space 



Hmut — Qrel ^^mat-Cmatf^mat Qrel with < 
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for the matter potential, it is precisely in this subspace where the weakest constraints can 
be placed. 

For the analysis of atmospheric, LBL and CHOOZ data, one can make the additional 
simplifying assumption of setting Am|^ = 0, hence neglecting the solar splitting. In this 
limit the 612 angle and the (5cp phase become unphysical, even in the presence of the gener- 
alized -ffmat introduced above. Thus altogether the relevant flavor transition probabilities 
for atmospheric, LBL and CHOOZ neutrinos depend on eight parameters: (Amg^^, ^13, 
^23) for the vacuum part, (e, 9912, (fiz) for the matter part, and (ai, a2) as relative phases. 

We notice that within this parametrization and neglecting Am^i both the vacuum and 
the matter part of the Hamiltonian have each two degenerate eigenvalues. In Ref. [53] it 
was noted that when this occurs the evolution equation can be reduced to an effective two- 
flavor problem in terms of an intermediate basis. Such effective 2x2 evolution Hamiltonian 
depends only on 3 parameters: Amg^^, e, and an angle which is a given function ^13, ^235 
(^12 and (/?i3 (see Appendix A for details). Therefore we can solve the evolution in this 
reduced parameter space and later on project back to the flavor basis. Technically this 
makes it possible to perform the analysis of the atmospheric, LBL and CHOOZ data in 
the full eight-dimensional parameter space. 

A theoretical framework for our proposed parametrization of the matter potential is 
provided by NSI affecting neutrino interactions in the Earth matter. They can be described 
by effective four-fermion operators of the form 

C^si = -2V2GFei^^{l7^j^'iyp)(f-ff,Pf) , (2.4) 

where / is a charged fermion, P = (L, R) and are dimensionless parameters encoding 
the deviation from standard interactions. Ordinary matter is composed by electrons (e), 
up-quarks (u) and down-quark (d), and in principle non-standard neutrino- matter inter- 
actions can involve any of these particles with different strength. In practice, however, 
the proton/neutron ratio found in matter is reasonably constant all over the Earth, and 
pretty close to 1. This implies that neutrino oscillations are only sensitive to the sum of 
these interactions, weighted with the relative abundance of each particle. Furthermore NSI 
enter in neutrino propagation only through the vector couplings = -|- e^^. We 
can therefore define the relevant combination of NSI parameters entering into the neutrino 
propagation in the Earth as: 

'-0^ E (^) 4/3 + 3£:i;3 + 3^/5 (2.5) 

f=e,u,d 

where Nf is the number density of the fermion /. The corresponding matter Hamiltonian 
reads: 

/I + £c.f, eer\ 
i?mat = V2GFNe{r) 6*^ 6^^ . (2.6) 

which includes both the standard (accounted by the term in the ee entry) and the 

non-standard interactions (accounted by the £q,^ terms) . 

The Hermitian matrix i^mat depends on 8 physical parameters. Oscillation experiments 
are sensitive to the matter potential up to an overall multiple of the identity, thus we can 
measure two differences of the £aa parameters (which must be real) as well as £e/i) Ser 
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and e^r (which, in general, can be complex). These 8 parameters are reduced to the 5 in 
Eq. (2.3) under the condition that two eigenvalues of ifmat are equal, which eliminates one 
real parameter, one angle and one phase leading to: 



£ee -^IJ.U= ^ (COS^ Lpi2 - SIT? 9912) COS^ (/3i3 - 1 , 



£ {sit? ipi^ - S\T? ipi2 COS^ (/Jis) , 

-£ cos(/3i2 sin 9912 cos^ (/9i3 e**-'^^"'^^^ 



(2.7) 



£er = ~^ COS (/9i2 COS 9913 Sm (fi^ € 



i{2ai+a2) 



£ sm ipi2 COS 9913 sm (^13 e 



i(ai+2a2) 



Notice that the term " — 1" at the end of Ege — accounts for the standard matter term. 
In terms of the NSI parameters e^^ the magnitude of £ can be always expressed as: 



£ = - mm 
2 \ 



+ £a0 — ^Saj: 



a/3 



. a0 a 



(2.8) 



and its sign can be defined by the relative size of the eigenvalues of i/mat^ £ is positive 
(negative) if the two equal eigenvalues are smaller (larger) than the third one. 

As mentioned above in Ref. [37] it was shown that strong cancellations (hence weak 
bounds on the e's) in the oscillation of atmospheric neutrinos occur if two of the NSI 
eigenvalues coincide. This is so because in this limit there is a two-dimensional subspace 
where i/^ — >• oscillations occur as in vacuum. Consequently, although Eq. (2.3) cover 
only a part of the most general parameter space for the matter potential induced by NSI, 
it is the part where the weakest constraints on NSI can be placed. Consequently our 
conclusions on the allowed values of the NSI parameters are not expected to be weakened 
once the "equal eigenvalues" approximation is removed. 

We conclude this section by noticing that until recently it was usually assumed that the 
bounds from loop contributions to flavor changing processes induced by the NSI implied 
£(,^ < 10"'^ [54] . If the condition = is imposed (together with the assumption that two 
NSI eigenvalues are equal) one is left with two possible scenarios: (i) £f^r = 0, which is a 
subset of the er sector studied in Refs. [37,38] (who made the assumption eg/i = g/iT = but 
left both eigenvalues free) , and (ii) Egr = 0, which corresponds to the NSI being dominantly 
in the /xr sector [3, 30, 36] and is revisited in Sec. 4. However, in Ref. [55] it was argued 
that the loop bound on £e^ does not hold in general and consequently it raised the interest 
of the general analysis here presented. 



3. Results 



Let us now present the results of a combined analysis of atmospheric, LBL and CHOOZ 
experiments in the context of oscillations with the generalized matter potential in the 
Earth in Eq. (2.3). It is well known that important constraints on flavor oscillations 
driven by Am|^ arises from the negative results on z/g disappearance at short baselines at 
the CHOOZ reactor experiment [56]. We include in our analysis the results of CHOOZ 
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under the assumption that, because of the short basehne, matter effects are irrelevant in 
this experiment. For atmospheric neutrinos we include the results of our reanalysis of the 
Super-Kamiokande I+II+III phases as presented in [57] and described in Ref. [9]; details 
on our simulation of the data samples and the statistical analysis can be found in the 
Appendix of Ref. [3]. For what concerns LBL accelerator experiments, we combine the 
results of the energy spectrum of events obtained by MINOS after an exposure to the 
Fermilab NuMI beam corresponding to a total of 7.25 x 10^*^ protons on target [58] with the 
data on antineutrino disappearance corresponding to 1.71 x 10^*^ protons on target [59]. 
We also account for the recent MINOS results on — )■ Ve transitions for which we fit the 
total event rate of observed events above the expected background based on an exposure 
of 7.01 X 10^° protons on target [60]. 

We present the results of the global analysis in Figs. 1, 2, 3 and 4. For simplicity the 
results are given for fixed values of the phases ai = a2 = 0. We have verified that the 
phases have little impact on the results of this analysis. In particular, no information on 
them can be extracted from the data {i.e., their allowed range is the full < < 27r 
interval) and they have only a minor effect on the determination of the other six real 
parameters. 

Fig. 1 contains the three two-dimensional projections of the allowed regions in the 
oscillation parameters jAmlj^j, sin^ 6*23 and sin^ ^13 after marginalizing over e, ipi2, and 
(^13. The regions are shown at 90%, 95%, 99% and 3cr CL (2 dof). Notice that since we 
have chosen to absorb the relative sign between the Am|;^ and the matter potential in e 
there is only one oscillation region. For the sake of comparison we also show in the figure 
the corresponding regions for the case of standard matter potential with normal ordering 
(e = +1, <^i2 = (^13 = 0) and inverted ordering {e = —1, (pi2 = (pi3 = 0). Fig. 2 displays 
the two-dimensional projections of the allowed regions in the matter potential parameters 
e, (pi2 and (/?i3 after marginalizing over lAmg^^l, ^23 and ^13. The regions are shown at 90%, 
95%, 99% and 3a CL (2 dof). From these figures we can draw the following conclusions on 
the determination of the matter potential: 

• No bound on the magnitude of the matter effects, e, can be derived from the analysis 
of atmospheric and LBL (and CHOOZ) experiments only in this general scenario. 
Consequently in order to attain bounds on e it is necessary to combine these results 
with those from other neutrino oscillation samples involving oscillations, such as 
solar and KamLAND experiments. However, a bound on e can be derived if a certain 
flavor structure of the matter potential is assumed a priori, implying that 9912 and/or 
</?i3 are larger than a certain amount. This is the case, for example, if we assume 
that no matter effects are present in the e/x and er projections, which corresponds to 
<^12 = 7r/2. We will study this case in the next section. 

• Matter potentials with |e| > 0(0.2) are only allowed as long as their flavor projections 
out of the ee entry are severely constrained. The constraint is stronger on rotations 
over the n flavor, (pu- We also notice that for large values of e the precise bounds 
on ipi2 and (pis present an oscillatory behavior. We will discuss the origin of this 
beaviour below. 

Qualitatively these results can be understood as follows. For ^/2GFNe{r)e 3> ^Tn\i/ {2Ei^) 
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Figure 1: Two-dimensional projections of the allowed regions from the global analysis of atmo- 
spheric, LBL and CHOOZ data in the oscillation parameters |Am§]^|, sin^ 6*23 and sin^ 6*13 after 
marginalizing over e, 1^912, and ipis (full regions). The best fit point is marked with a star. For the 
sake of comparison we also show in the figure the corresponding regions for the case of standard 
matter potential with normal ordering (e = +1, ipi2 = fis — 0) as unfilled regions with full lines 
with best fit point marked with a circle and inverted ordering (e — —1, ipi2 — ipis — 0) as unfilled 
regions with dashed lines with best fit point marked with a triangle. The regions are shown at 90%, 
95%, 99% and 3cr CL (2 dof). 

the neutrinos can undergo flavor oscillations with two different phases: 

Avac = X /(023, ^13, <^12, (3.1) 

^ _ V2GFN,ir)eL 

^mat — 2 l-J-^J 

where /(023) ^13; ^12, fis) is a combination of trigonometric functions of the four angles: for 
example, for ^13 = we get / = cos^ ipi2 cos^ ipis + (cos (^13 sin (^12 cos ^23 — sin 9913 sin023)^- 
Ufj^ o transitions oscillate dominantly with phase Avac^ this is the generalization of the 
results in Ref. [37,38] and it means that in this scenario — >• Ur oscillations have the same 



-7- 




e 



Figure 2: Two-dimensional projection of the of the allowed regions from the global analysis of 
atmospheric, LBL and CHOOZ data in the matter potential parameters e, 1^912, and (^13 after 
marginalization with respect to the undisplayed parameters. The regions are shown at 90%, 95%, 
99% and 3cr CL (2 dof). The best fit point is marked with a star. 

dependence on the neutrino energy and distance as vacuum oscillations, even for large e, 
which opens the possibility of giving a good description of atmospheric and LBL results. 
In the same regime o fe and -f-)- z^r transitions proceed dominantly via oscillations 
of phase Amat 

Pefj. = cos^ (/?i3 sin^(2(/7i2) sin^ Amat , (3.3) 
Per = cos^ ipi2 sin^(2((C)i3) sin^ Amat • (3.4) 

This means that even for ^13 = atmospheric fg's can disappear into either or Ur 
and they do so independently of their energy. Similarly, some of the atmospheric muon 
neutrinos will oscillate into electron neutrinos independently of their energy. This is in 
clear conflict with the atmospheric data and therefore the amplitudes of these oscillations 
are constrained. This leads to the strong bounds on |(/Ji2| and {(fisl observed in Fig. 2 for 
|e| > 0(0.2). Furthermore, since the atmospheric Vg data mostly constrain the sum of 
Pe^j. + Per, the bounds on ipi2 and ipis for a given e are strongly correlated. Notice also 
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that atmospheric neutrino angular bins correspond to given ranges of L, rather than just 
a specific value as in the case of MINOS, so for large e the phase is averaged in L and 
the atmospheric bounds on (pi2 and ipis become independent of e. 

Similarly to the atmospheric neutrino data, the MINOS t'^ disappearance spectrum 
is well explained by the vacuum-like oscillations — t- Ur- But it also leads to further 
constraints on the contribution from Pg/* in (3.3), since the presence of this term would 
induce an energy independent contribution to the survival probability which is in conflict 
with the observations. This tightens the atmospheric constraint on ipi2- Since L in MINOS 
is fixed, the contribution of Pg^i to the MINOS spectrum has an oscillatory behavior on 
e (i.e., the phase is not averaged over L as long as the size of the detector is small 
compared to [GFNe{r)£]~^) which leads to the "spike-like" shape of the allowed region 
for (pi2 versus e. Given the strong correlation between the allowed values of cpi^ and 9912 
imposed by the atmospheric data, this "spike-like" shape is also projected in the allowed 
region for (^13 versus e. 

For what concerns the standard oscillation parameters, our analysis shows that their 
present determination is rather robust even in the presence of a generalized matter potential 
of the form (2.3). To further illustrate this we show in Fig. 3 the allowed regions of the 
three oscillation parameters as a function of the matter potential magnitude e. As seen 
in the figure, the allowed ranges for the oscillation parameters do not get substantially 
modified for any value of e (see also Eq. (5.1)). 

The dependence on these results on the data samples included in the analysis is shown 
in Fig. 4 where we plot the one dimensional projections of the Ax^ for three different 
combination of experiments as a function of each of the six parameters. In each case 
the Ax^ is defined with respect to the minimum obtained for the standard oscillation case 
which occurs in the normal ordering (e = +1, (pi2 = V'la = 0). We read from the figure that 
for any of the three combination of data shown, there is no substantial improvement of the 
fit due to the inclusion of the extra parameters. Ax^ is reduced by a maximum of 1.5 units 
in the global analysis. In other words, even though no constraint on the magnitude of the 
matter effects can be placed by the data, the data do not show any need of a non-standard 
matter potential. In particular comparing the results including and not including the z?^ 
MINOS results we see that inclusion of non-standard matter effects does not result into a 
substantial improvement on the quality of the combined fit. Prom the lower panel of the 
figure we sec that inclusion of the MINOS antincutrino results implies the appearance of a 
new local minimum at e = —3. This local minimum occurs at (pi3 = —16° and (pi2 = —7° 
and it is the equivalent of the best solution found in the partial analysis of Refs. [61-63] 
to explain the discrepancy of MINOS neutrino and antineutrino results in the framework 
of NSI. However, as seen in the figure, we find that the "improvement" on the quality of 
the description associated with this new minimum is about 1 unit in the global x^- This is 
partly due to the constraints from atmospheric neutrino data on the required values of the 
NSI and partly on details associated with the statistical analysis of MINOS in the presence 
of NSI, which we will quantify in more detail in the next section. 

As a matter of curiosity we notice that the best fit of the global analysis corresponds 
to "standard" magnitude matter effects with inverted mass ordering (e = —1), maximal 
023 and 9i3 = 0, but with non-standard flavor projections (^12 = 8° and (^13 = 20°. This 
behavior is already present when only atmospheric and CHOOZ data are included in the 
analysis and it is driven by the small excess of i/^ events in the multi-GcV atmospheric 
samples. This excess was more significant in SK-I and was highlighted in Ref. [64] as a 
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Figure 3: Two-dimensional projection of the of the allowed regions from the global analysis of 
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with respect to all other parameters. The Ax^ is defined with respect to the minimum for pure 
oscillations which lies in the e = +1, </?i2 = ifiis = subspace. 
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possible "hint" for non-zero ^13, although as pointed out in Refs. [9, 11] the statistical 
significance of this effect was reduced in SK-II and SK-III samples. What we find in our 
analysis is that it is slightly easier to fit the data with a non-standard flavor structure of the 
matter potential and a zero ^13 than with standard matter potential and a non-vanishing 
^13. Nevertheless the difference, yielding a Ax^ = —1.5, is not statistically significant. 

Finally we show in Fig. 5 the two-dimensional allowed regions projected in terms of 
the equivalent NSI parameters SajS, see Eq. (2.7). The figure illustrates that although no 
absolute bound can be place on the individual Eq,^ parameters, they can only be large if they 
are highly correlated among them. In other words, in order to comply with the observed 
data in atmospheric, LBL and CHOOZ experiments, a beyond the SM scenario which could 
generate such large NSI parameters must have a very constrained flavor structure which 
implies, for example, Sf^r/s ^ 0.1 which, unfortunately, is not the one required to explain 
the MINOS Un versus i>ij, results. We elaborate more on this in the following section. 

4. No matter effects in e/x and er sectors 

In this section we study the particular case in which no potential is present in the e// and 
er sector, or, equivalently, that the NSI parameters eg/* = £er = are known to be severely 
constrained from other sectors. As described in Refs. [30,36], in this case a bound on the 
strength of the potential arises from the analysis of atmospheric and LBL neutrinos. From 
Eq. (2.7) we see that Cg^ = £er = corresponds to (pi2 = tt/2. If we further assume ^13 = 
then the Ve decouples from the evolution equation. Moreover one of the two phases become 
unphysical so that the matter part of the Hamiltonian depends on three parameters. In 
Refs. [3, 36] we labeled them T, ^ and rj and they are related to the present notation by 

J'=]Im+^- ''~/^^^' = '^, e=^=</'i3, r? = arg(£^,) = ai + 2a2 . (4.1) 
The relevant oscillation probabilities are given by 

(Am? L \ 

where 

sin^ 2G = ^ (sin^ 26'23 + sin^ 2(^i3 -|- 2R sin 26*23 sin 2(pis cos 77) , 
TZ = + + 2R (cos 2^23 cos 2(^i3 -|- sin 2^23 sin 2(/?i3 cos rj) , 

2 A^^ii 

So in this case, for large e the unique oscillation wavelength becomes energy independent 
in contradiction with the data, and an upper bound on the magnitude of the non-standard 
matter effects can be placed. Also, since in this case we are fixing (^12 = 7r/2 the standard 
two-neutrino — )• Vr oscillation scenario is recovered in the limit e = 0. 

In Fig. 6 we show the two-dimensional projections of the allowed parameter region 
for the analysis of atmospheric and LBL data. The filled areas in the upper left panel 
show the allowed regions on the (Am|]^, sin^023) plane. For the sake of comparison we 
also show the lines corresponding to the contours in the absence of new physics. We see 
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Figure 5: Two-dimensional projection of the of the allowed regions from the global analysis of 
atmospheric, LBL and CHOOZ in terms of the NSI parameters Sap versus the e parameter. The 
regions are shown at 90%, 95%, 99% and 3a CL (2 dof). The best fit point is marked with a star. 
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Figure 6: Allowed parameter regions for the analysis of atmospheric and LBL data in presence 
of — > Vr oscillations and non-standard matter effects. Each panel shows a two-dimensional 
projection of the allowed five-dimensional region after marginalization with respect to the three 
undisplayed parameters. The different contours correspond to the two-dimensional allowed regions 
at 90%, 95%, 99% and 3a CL (2 dof). The filled areas in the upper left panel show the projected 
two-dimensional allowed regions on the {Amfi, sin^023) plane. The best fit point is marked with 
a star. For the sake of comparison we also show the lines corresponding to the contours in the 
absence of new physics and mark with a triangle the position of the best fit point. The regions in 
other panels show different projections of the allowed values of the parameters characterizing the 
strength and mixing of the new matter effects. The full regions correspond to the case r] S {0, tt} 
with best fit point marked by a star, while the solid lines correspond to arbitrary values of the 
phase 77 with best fit point marked by a triangle. In the two lower panels we show for comparison 
the corresponding regions obtained in the general analysis in Fig. 5 (dashed contours). 
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that in this constrained scenario the determination of the oscillation parameters is very 
much unaffected by the presence of the non-standard matter effects. The other panels of 
Fig. 6 show different projections of the allowed values of the parameters characterizing the 
strength and mixing of the non-standard matter effects. The solid contours correspond to 
arbitrary values of the phase rj, while the colored regions correspond to the case rj € {0, tt}. 
Translated in terms of NSI the corresponding 90% {3a) bounds read: 

(-0.055) -0.035 < £ij,r < +0.018 (-F0.035) , \£rr - ^ 0-097 (0-16) (4.6) 

for real NSI, and 

k/xrl < 0.035 (0.055) , \err - £mmI ^ 0.11 (0.18) (4.7) 

for the general case of complex e^,-. For the sake of comparison we show in the two lower 
panels of Fig. 6 the corresponding regions obtained in the general analysis in Fig. 5 as 
dashed contours. We see how only in the restricted analysis a bound on the e parameter is 
set. However the corresponding flavor projections are not more severely constrained in this 
restricted analysis when compared to the general one presented in the previous section. 

To finish this section we turn to the possible role of non- vanishing e^r and Ett — ^nfj. 
to resolve the tension between the oscillation results for neutrinos and antineutrinos in the 
MINOS experiment [58,59]. As shown in the previous section, even in the most general 
case we found that the inclusion of the matter potential leads to a decrease of the Ax^ of 
at most ~ 1.5 units with respect to pure oscillations. To further quantify the effect of NSI 
in the fir sector in the MINOS data we show in Fig. 7 (right panel) the Ax^ of the analysis 
of MINOS neutrino and antineutrino results in the framework of V/^ Vr oscillations in 
the presence of non-vanishing NSI in the jur sector. For the sake of comparison with the 
literature we show the results for two analysis, one in which all the data of all energy bins 
of the MINOS neutrino and antineutrino spectrum are included and one where only the 
bins with 1 < E',^ < 5 GeV and 1 < E'p < 8 GeV are included as in Ref. [61]. We also show 
in the left panel of Fig. 7 the corresponding oscillation regions for the analysis of MINOS 
in the framework of pure Am^i — )■ Ut oscillations for the two analysis. As seen in the 
figure, in the framework of Am3^ oscillations only both analysis yield equivalent oscillation 
regions. This is so because for pure Am\i oscillations the data of the higher energy bins has 
little relevance in the determination of the oscillation parameters since at those energies 
the Aml^-driven oscillation wavelength is much longer than the characteristic L in MINOS 
so the corresponding survival probability is 1. However, when NSI are included, the and 
Vfj, survival probabilities acquire an energy independent piece as seen in Eq. (4.2). As a 
consequence, in the presence of NSI the inclusion of the high energy bins is relevant to the 
conclusions of the analysis. This is quantified in the right panel of Fig. 7 where we show the 
^^MINOS five-parameter analysis (in terms of oscillations plus non-standard matter 

effects in the jiT sector) as a function of the strength parameter for the non-standard matter 
potential after marginalizing with respect to the other four parameters (Am^^, ^23, Viz-, v)- 
Here Ax^jj^Qg is defined with respect to the pure oscillation scenario. From this figure we 
read than when only the low energy bins are included, NSI seem to be able to improve the 
agreement between and i/^ disappearance data by Ax^^jj^Qg ~ —7, in good agreement 
with the results in Ref. [61]. However the inclusion of the effect of these same NSI in the 
high energy bins spoils this improvement and leads only to Ax^jj^os — Furthermore 
the required strength for the NSI parameters is in conflict with the atmospheric results, as 
shown in the previous section and correctly noted in [61]. 
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Figure 7: Results of the analysis of MINOS i'^ and Df^ disappearance data for two sets of analysis 
differing in the inclusion of the high energy bins (see text for details). The left panel show the 
allowed regions in the framework of pure Am^j^-driven oscillations. The full regions are the two- 
dimensional allowed regions at 90%, 95%, 99% and Scr CL (2 dof) for the analysis including all the 
bins and the corresponding best fit point is marked with a star. The unfilled regions correspond 
to the analysis where only the 1 < < 5 GeV and 1 < i?p < 8 GeV bins are considered (best fit 
point marked with a triangle). The right panel show the Axminos ^ ^ function of the NSI strength 
for the five-parameter analysis (in terms of oscillations plus non-standard matter effects in the fir 
sector) after marginalizing with respect to to the four undisplayed parameters. The lower (upper) 
curve corresponds to the analysis without (with) inclusion of the high energy bins. 



5. Summary 

In this work we have addressed the question: what do atmospheric and LBL neutrino data 
teach us about the presence and flavor structure of matter effects in propagation? 

First, to set out to answer this question we have introduced in Sec. 2 a phenomeno- 
logical parametrization of the matter potential allowing for departures of the SM one by 
rescaling of its strength, rotation from the ee sector, and rephasing with respect to vacuum. 
In this scenario after setting Am^]^ = in the analysis of atmospheric, LBL and CHOOZ 
data, the relevant flavor transition probabilities depend on eight parameters: the three 
oscillation parameters (Amgj^, 623 and ^13), three matter potential parameters (e which 
represents a rescaling of the matter potential strength, (^12 which allows for projection of 
the potential into the flavor, and ipis which allows for its projection into the I'r flavor) 
and two relative phases. The relation to the matter potential induced by NSI is given in 
Eq. (2.7). 

The answer to the question in the most general case is given in Sec. 3. We flnd that 
in general the strength of the matter potential e cannot be determined by the analysis. 
This is the generalization of the results in Ref. [37, 38] and it is due to the fact that for 
\/2GFNf,{r)e ^ Am|]^/(2£'i,), o transitions oscillate dominantly with the phase 
Avac in Eq. (3.1) and have the same dependence on the neutrino energy and distance as 
vacuum oscillations. This allows for good description of atmospheric and LBL results even 
for very large values of e. 
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However important information is obtained from the data analysis on the flavor struc- 
ture of the matter potential. This is so because for large e the -H- (Eq. (3.3)) and 
-H- Vt (Eq. (3.4)) transitions proceed dominantly via oscillations with the energy inde- 
pendent phase Aniat in Eqs. (3.2) and amplitudes controlled by the projection angles ipi2 
and 9913. This means that even for ^13 = atmospheric z^g's can disappear into either 
or i/t and they do so independently of their energy. Similarly, some atmospheric and 
LBL f^'s will oscillate into electron neutrinos independently of their energy. This is in 
clear conflict with the atmospheric and LBL data, and as a consequence potentials with 
|e| > 0(0.2) are only allowed as long as their flavor projections out of the ee entry are 
severely constrained. As can be seen in Fig. 2, the constraints are stronger on rotations 
over the jx flavor, controlled by ip\2- 

We also find that the present determination of the oscillation parameters is rather 
robust even in the presence of this general form of the matter potential as illustrated in 
Figs. 1 and 3. For example the 3a" ranges of the oscillation parameters read: 

Oscillations Oscillations + Generalized Matter Potential 

2.11 X 10-^ < |Am|i| < 2.69 X 10-3 2.10 x 10"^ < |Amii| < 3.06 x 10"^ 

0.34 < sin2 ^23 < 0.68 0.26 < sin^ ^23 < 0.68 ^ ' ' 

sin^ < 0.051 sin^ < 0.060 

Concerning the quality of the fit, we find that even though no constraint on the mag- 
nitude of the matter effects can be placed by the data, the data does not show any favoring 
for a non-standard matter potential. There is no substantial improvement on the statistical 
quality of the fit due to the inclusion of the extra matter parameters with A^^ reduced by 
a maximum of 1.5 units in the global analysis. This is so even after the inclusion of the 
MINOS disappearance results. Thus we conclude that within the context of a global 
analysis the inclusion of the NSI does not significantly alleviate the tension between the os- 
cillation results for neutrinos and antineutrinos arising from the MINOS disappearance 
results. This is partly due to the constraints from atmospheric neutrino data and partly 
to the fact that the non-standard matter potential parameters required to better fit the 
MINOS energy spectrum of and events with < 5 GeV and Ey <S GeV with the 
same oscillation parameters worsens the description of the spectrum of Ei, > 5 GeV events 
as shown in the right panel of Fig. 7. 

Finally we have revisited the particular scenario in which matter effects in the ejj, and 
er sectors are very suppressed (Sec. 4) . In this case a bound on the strength of the NSI in 
the /ir sector can be set. The corresponding updated bounds shown in Eqs. (4.6) and (4.7) 
are 5-20% stronger than our last published bounds in Ref. [3]. 
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A. Parametrization of the matter potential 



We give here some details on the simphfications associated with the condition that both 

the vacuum and the matter part of the Hamiltonian have each two degenerate eigenvalues. 
In the most general case the matter potential can be parametrized as: 



^^^mat = QrelC^mat-DmatC^iat^Ll "^^^^ ' 



' Q,el = diag(e^°^e^"^e-^"l-^"2), 
I L»mat = V2GFiVe(r) diag(e, e', 0) . 



Setting e' = implies that the ip^z angle and the 5ns phase become unphysical. In the 
limit Am|i = the general Hamiltonian reduces to: 



n — n23-n,]^3i>'vac -0-13-0-23 "r Vrel-f'-i2-0'i3-'->'mat-n-i3-n.]^2 Vrel 

where Rf, = Rij{6ij) and Rf- = Rij{(pij). We can write: 



(A.2) 



H — -ft23-Ki3 



Dvac + ^-Rl3-R23Qrel-Rr2-Rr3^ -^mat ^-^13 -^12 '5i[el-^23-^13) -^13-^23 i-^-^) 



where we simply collected the two real rotations i?23 and R^^. Changing the parametriza- 
tion: 

■^13-^23Qrel-?^12-^13 ^ QeflF-f^i2-^l3-^23 Qmaj (-^-4) 

where Qeff = diag (e'^^, e*''^^ g-«/3i-«/32^ relative phases, i?f2 and i^fg are real 

rotations by angles ^12 and ■^13, i?!/ is a complex rotation by angle '023 and phase 623, and 
Qmaj contains two Majorana phases. The matrix on the l.h.s. is a det = 1 unitary matrix, 
so it can be decomposed as in the r.h.s.. The new angles (V'i25 V'lS) ^23) as well as the 
phases (523; Pi, h) and the Majorana phases in Qmaj are some complicated functions of 
the original angles (^13, ^23) ^12-, V^is) and of the phases (cci, 02), but they do not depend 
on Am|^, (neutrino energy), e, or Ne{r). They are universal and independent of the 
details of the trajectory. 

Now, obviously the Majorana phases Qmaj disappear, as well as 623 and ^23 (since they 
rotate the (0, 0) block in Draat)- Hence we are left with: 



H — -tt23-n-13 

Now we note that: 



Dvac + Qeff-Ri2-^l3-^mat-?^i3 -^12 Qeff 



Dyac — QefF-f^i2-^vac-^i2 Qcff 



-^13-^23 



-'--'vac "r -n,]^3J-y'niat-''-i3 



where U,e = R^23Ri3QesRt2 



since i?^2 rotates the (0, 0) block in -Dyac and QeS is diagonal. Therefore: 

H = U,s 
This lead to a scattering matrix: 

S = UesSesUlg with S'eff = evolution-of ^Z)vac + -Ria-Dmat-Ril) 



(A.5) 
(A.6) 

(A.7) 
(A.8) 
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It is immediate to see that the effective Hamiltonian H^s = ^vac + ^la-Dmat^iJ depends 
only on 3 parameters (Am^i, e, tpis) and has the form (up to an overaU multiple of the 
identity): 

/Hd Hn \ (Sd -S*^\ 

i^eff = i^p =^ 5eff = 5p (A.9) 

\Hm -HdJ \Sn S}, I 

hence it factorizes into a 2 x 2 block and a 1 x 1 block. The remaining 5 parameters ^13, 
^23) /?2, V'12 are non-dynamical and can be reinserted later to form the x^. Concretely, 
the neutrino probabilities can be written as the sum of 18 different terms: 

18 

P.^ = Y.CfPf (A.IO) 

n=l 

where = cf {613,623, Pi, P2,M and Pf = Pf{Amli,e, V'ls) are real numbers re- 
lated to UeS and SeS by the following formulas: 





cf = 




pf = 


1, 




cf = 




peff _ 




cf 




wf*wf - 




pf + ipf = 


2 Su S*j\j , 


'-'5 


+ ^cf = 






pf + ipf = 


2 S"/) S*p , 


cf 


+ iCf = 






pf + ipf = 


2Sn Sp , 


^9 


+ ^Cf = 








2Sd Sd , 


^11 


+ ^Cti = 






pcff 1 „• pcff _ 


2Sn Sn , 




+ iCf = 


w^fw^i - 


wf*wf, 


peff 1 „• pcff _ 
^13 + ^^li — 


2Sd Sn , 


^15 


+ ^< = 






peff 4_ peff _ 


2 iSp , 










peff 4_ peff _ 


2Sn Sp , 



(A.ll) 



with Wjj- = U^j . For antineutrinos one has to flip the sign of the matter term in ifeff 
and replace UeS U*q (or, equivalently, SeS — >■ S*g) in Eq. (A.ll). 



B. Symmetries and pairameter ranges 

In order to describe the physical ranges for the parameters in the problem we study which 
kind of transformations leave the probabilities unaffected. Let's consider a problem de- 
scribed by an Hamiltonian H. The transformation H — t- —H* implies exp{—iHL) — )■ 
exp{iH*L) = exp[{—iHL)*] = [exp(—iHL)]*, hence it leaves the probabilities unchanged. 
Such transformation is realized by: 

(0) Amli —Am^i A e ^ —e A ai ^ —ai A 02 ^ —02; 

This symmetry implies that only the relative sign of Am^i and e matters, so that we can 
either assume e > and study both signs of Arn^i, or assume Am^i > and study both 
sign of e. Furthermore, any rephasing H — )• QHQ* where Q = diag (e*", e*'', e**^) leads 
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to a rephasing of the scattering matrix exp{—iHL) — )■ Q exp{—iH L)Q* , which leaves the 
probabihties unaffected. Our Hamiltonian can be written as: 

H = U^^cD^^cUL + Qrelf/matl^matf/LtQll (B-l) 

where [/vac = ^23(^23)-Ri3(6'i3) and t/mat = Ri2{(pi2)Ri3{(pi3)- A transformation H 
QHQ* is achieved if: 

Uvsic QUvacQ' and Qrcl^^mat QQidUmatQ" (B.2) 

wlicrc Q, Q' , Q" are all pure phase diagonal matrices. Hence, any transformation of 
the parameters which is equivalent to a specific choice of Q, Q' and Q" will leave the 
probabilities unaffected. With this, we can generate a list of symmetries: 

(1) ^ TT + ^13. This is equivalent to Q' = diag(— 1, +1, —1); 

(2) 6*13 -(9i3 A 623 -^7T + 023- This is equivalent to Q' = diag(+l, -1, -1); 

(3) 9913 — )• TT + ipi3. This is equivalent to Q" = diag(— 1, +1, —1); 

(4) ifis —(pi3 A ipi2 TT + ipv2- This is equivalent to Q" = diag(— 1, —1, +1); 

(5) 023 -023 A ipi2 -(^12- This is equivalent to Q = Q' = Q" = diag(+l, -1, +1); 

(6) 023 ^ 7T — 023 A ipi2 ^ TT + ipi2. This is equivalent to Q = diag(+l, +1, — 1), 
Q' = diag(+1,-1,+1), Q" = -I; 

(7) ipi2 — > —y^i2 A ai -> ai + 7r/3 A 0:2 — >■ 0:2 — 27r/3. This is equivalent to Q" = 
gW3 diag(+l,-l,+l). 

(8) (pi2 7r+(^i2 A ai ^ ai+Tv/S A 02 ^ a2+7r/3. This is equivalent to Q" = e~^'^/^/. 

With these symmetries wc can reduce the range of the parameters (^23, ^13, fi2, fis) from 
the most general [0 27r] to something less: 

• using (1) we can reduce ^13 to the range [— 7r/2 +7r/2]; 

• using (2) we can reduce ^13 to the range [0 — >■ 7r/2]; 

• using (5) we can reduce ^23 to the range [0 — > tt] ; 

• using (6) we can reduce ^23 to the range [0 — >■ 7r/2]; 

• using (3) we can reduce (^13 to the range [— 7r/2 Tr/2]. 
For the case of real NSI: 

• using (4) we can reduce (^12 to the range [—it/ 2 — t- vr/2]; 

so that for ai = ^2 = we can reduce the parameter range for the angles to < 0ij < 7r/2 
and — 7r/2 < tpij < tt/I. For the case of complex NSI: 

• using (4) we can reduce 9913 to the range [0 — > vr/2]; 

• using (8) we can reduce 1^12 to the range [—it/ 2 — )• 7r/2]; 

• using (7) we can reduce (pi2 to the range [0 7r/2]; 

so that in the general case of unconstrained ccj we can reduce the parameter range for the 
angles to < 0ij < 7r/2 and < ipij < 7r/2. 
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